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Abstract. - The ageing of the bond-disordered two-dimensional Ising model quenched to 
below its critical point is studied through the two-time autocorrelator and thermoremanent 
magnetization (trm). The corresponding ageing exponents are determined. The form of the 
scaling function of the trm is well described by the theory of local scale-invariance. 



A spin system quenched to below or exactly at its critical temperature T c undergoes 
ageing, that is time-translation invariance is broken and observables such as correlation or 
response functions display dynamical scaling. While systematic studies of ageing appear to 
have been carried out first in studies of the mechanical response of glassy systems [1], it has 
been realized relatively recently that a similar phenomenology also applies for simple ferro- 
magnets without disorder, see [2-5] for reviews. In these simpler systems, for both T = T c 
and T < T c there is a single time-dependent length-scale growing with time as L(t) ~ t x l z 
where z is the dynamical exponent. L(t) is identified with the typical linear size of corre- 
lated or ordered clusters, respectively. It is common to study ageing through the two-time 
autocorrelation function C(t,s) = (<f>(t,r)<f)(s,r)) or the two-time linear autoresponse func- 
tion R(t,s) — S(4>(t,r))/Sh(s,r)\f l=0 . Here <j>{t,r) denotes the space-time-dependent order- 
parameter, h(s, r) is the space-time-dependent conjugate magnetic field, t is referred to as 
observation time and s as waiting time. For ageing systems, averages depend on both t and s 
and not merely on the difference r = t—s. In the ageing regime i,s> i m icro and t— s 3> t m i cr o, 
where t m i cro is some microscopic time scale, dynamical scaling holds (at equilibrium a = b) 

C(t, s) = s- b f c (t/s) , R(t, s) = s- l - a f R (jt/s) (1) 

such that the scaling functions fc,n(y) satisfy the following asymptotic behaviour 

fc(y) - y- Xc/z , f R (y) - 2T Wz (2) 
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as y — y oo and where Ac and A/j, respectively, are known as the autocorrelation [6,7] and 
autoresponse exponents [8]. For non-conserved dynamics (model A), these exponents are 
independent of the equilibrium exponents and of z [9]. In simple magnets without disorder 
and with short-ranged initial conditions one finds Ac = A,r, see [2, 10] and this has been 
reconfirmed in a recent second-order perturbative analysis of the time-dependent Ginzburg- 
Landau equation [11]. On the other hand, for either long-ranged initial correlations in ageing 
ferromagnets [8] or else for disordered systems such as the random-phase sine-Gordon model 
[12], Ac and A# are known to be different from each other. Finally, if T < T c , then for simple 
magnets with short-ranged equilibrium correlation functions (such as the Ising model in d > 1 
dimensions) one has 6 = and standard scaling arguments show that a = 1/z. On the other 
hand, if T = T c , then a = b = (d — 2 + r))/z where 77 is a standard equilibrium critical exponent. 

Furthermore, according to the theory of local scale-invariance (LSI) [13] the response 
functions of ageing systems with an algebraic growth law L{t) ~ t 1 ^ are expected to transform 
covariantly under local space-time transformations t 1— ► (1 + s) z t, r t—y (1 + e)r with an 
infinitesimal e = s(t, r). Restricting the admitted transformations such that time-translations 
are excluded (as is natural to do in the context of ageing) leads to the following prediction of 
the scaling function fn(y), up to normalization [10,13-15] 

fR(y) = y 1+a '- XR/z (y-iV 1 - a ' (3) 

where a' is a new exponent. Over the past few years, evidence confirming the prediction @ 
has been accumulating. First, consider phase-ordering kinetics with a non-conserved order- 
parameter (where z = 2 [16]). Then the prediction @ with a — a! was confirmed in the exactly 
solvable random walk and the spherical model [17] and numerically in the 2D/3D Ising [18] 
and XY [19] models and the 2D g-states Potts model with q = 2, 3,8 [20]. Second, consider 
non-equilibrium critical dynamics with a non-conserved order-parameter, where in general 
z 2. Again, in many spin systems eq. © is either exactly reproduced in solvable models 
or describes very well the numerical data in in a large variety of systems, in general with a 
and a' being distinct, including ferromagnets such as the Glauber- Ising model in d = 1,2,3 
dimensions, the OJK approximation, the spherical model and glassy systems such as the 
Frederikson- Anderson model and the critical Ising spin glass, see the references quoted in [15]. 
In addition, ferromagnets in restricted gometries have also been studied [21,22]. Finally, tests 
of eq. (J3J) for ID critical systems without detailed balance (where holds for connected 
correlators with a ^ b in general, see [23-26]) were performed for the contact process (for 
t/s > 1.1) [15,23,27], the non-equilibrium kinetic Ising model [26] and bosonic models [28]. 
We stress that both the theoretical prediction © of LSI as well as almost all existing tests of 
it( 1 ) assume an initial state with non- vanishing magnetization. Finally, stochastic processes 
such as the zero-range process [29] lack the spatial structure assumed in local scale-invariance 
(for instance Galilei-invariance for z = 2 or a convenient generalization for z ^ 2) and do not 
reproduce 

In this letter, we study the ageing properties of the 2D Ising model with bond disorder, 
but without frustrations, which describes the ageing of diluted magnets. The hamiltonian is 

H = - ^ ■l,,n,G , ; cr.j = ±1 (4) 

where the random variables are equally distributed over the interval [1 — e/2,1 + e/2] 
with < e < 2. We consider a square lattice with 300 x 300 spins, which is large enough 

( )In the ID critical contact process a full lattice was used as initial state. It is conceivable that this 
distinct initial condition might cause the departure seen from eq. in the region t/s < 1.1 in the contact 
process [15,27] but the question is still open. 
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Table I - Quantities describing the ageing behaviour of the disordered 2D Ising model. 



e T 


b a ^c/z Xr./z 


r rx 


1/z 


0.5 1.0 

a q 
U.o 

0.6 
0.4 


0.005(2) 0.398(5) 0.570(5) 0.61(1) 

a aai a qqo /i ^ r\ kkc\( k\ acak/'ia\ 
U.VUi(Z) U.oo/^4J U.OOLHoJ U.oyo(lUj 

-0.008(2) 0.353(4) 0.520(5) 0.58(1) 
-0.010(2) 0.310(5) 0.460(5) 0.52(1) 


0.021(1) 

0.022(1) 
0.029(2) 0.008(1) 


0.400 
U.ool 
0.353 
0.308 


1.0 1.0 
0.8 
0.6 
0.4 


-0.014(2) 0.330(5) 0.490(5) 0.51(1) 
-0.014(2) 0.308(4) 0.450(5) 0.49(1) 
-0.015(2) 0.277(6) 0.380(5) 0.46(1) 
-0.013(2) 0.22(1) 0.290(5) 0.375(10) 


0.021(1) -0.009(1) 
0.019(1) -0.006(1) 
0.020(1) -0.007(1) 
0.026(2) -0.014(3) 


0.333 
0.308 
0.273 
0.222 


2.0 1.0 
0.8 
0.6 


-0.021(1) 0.24(2) 0.320(5) 0.33(1) 
-0.018(1) 0.22(2) 0.270(5) 0.30(1) 
-0.012(1) 0.17(2) 0.220(5) 0.27(1) 


0.048(2) -0.048(4) 
0.093(3) -0.042(4) 
0.194(4) -0.033(3) 


0.250 
0.222 
0.188 



such that no finite-size effects are notable. It is known that this model has a second-order 
phase-transition and that the value of of the critical temperature T c = T c (e) does not change 
much with e, since (Jij) = 1, such that 2 < T c (e) < 2.269 [30]. The kinetics (non-conserved 
order-parameter) is described through a standard heat-bath algorithm. Throughout, we use a 
fully disordered initial state and average over typically 10000 different realizations of the bond 
disorder and of the thermal noise. In this work, we shall limit ourselves to quenches deep into 
the ordered phase T < T c (e). Recently detailed studies were carried out by Paul, Puri and 
Rieger [30,31] on the growth law in the low-temperature phase. Assuming that the energy 
barriers for domain- wall motion depend logarithmically on the linear domain-size L = L(t), 
they proposed for a non-conserved order-parameter the growth law, valid for t — > oo [30] 

L(f)~W 2 , z = z(T,e) = 2 + e/T (5) 

and they confirmed this through extensive simulations of the single-time correlation functions. 
It will be one of our objectives to further test this relationship. In addition, we see that the 
value of the dynamical exponent can be easily changed so that tests of the LSI-prediction 

for a large range of values of z becomes possible. This is particularly interesting since in 
previous tests of LSI in nonequilibrium critical dynamics the value of z always remained quite 
close to 2, the largest deviations occurring in the ID contact process, where z ~ 1.6. 

We begin the discussion of our results by considering the two-times autocorrelation function 
C(t,s). In figure^, we plot our data taken for e = 2 and T = 1 of the autocorrelator 
over against t/s. In constrast to what is found in phase-ordering for simple, non-disordered 
magnets, no data-collapse occurs here. Comparing with the anticipated long-time scaling 
lfT)l. this implies that the exponent b cannot vanish, if that scaling form is valid at all. In 
figure ^) we show that a relatively good, but not perfect scaling can be achieved with a non- 
vanishing effective value of b. Still, in spite of the large waiting times used, apparently there 
remain important finite-time corrections to the anticipated scaling form ifTjl. A forthcoming 
study [32] of the site-disordered Ising model rather suggests C{t, s) ~ s~ b (*/ s ) fc(t/ 's) with two 
continuous functions fc and b. Phenomenologically, this certainly achieves a better collapse 
than with a constant b. More detailed studies on the scaling of C(t, s) will be needed. 

Assuming the standard scaling form (11121) being valid for C, we extract the value of the 
autocorrelation exponent Xc/z, listed together with effective values for b in table ||| 

We now turn to a discussion of the scaling of the linear response. While a direct calculation 
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Fig. 1 - Autocorrelation for e = 2 and T — 1 as a function of t/s for various waiting times s. For 
b = no data collapse is observed (a), whereas for b = —0.021 the data collapse nicely for large values 
ott/s (b). 



of the functional derivative for i? would produce extremely noisy data, it is of interest to 
consider an integrated reponse function [33] where much of the noise is smoothed out by 
the integration. Here we work with the thermoremanent magnetization Mtrm(A s) which is 
obtained after applying a small spatially random magnetic field with amplitude h = 0.05. 
From our experience with the scaling of Mtrm(£, s) in the phase-ordering of simple magnets, 
we expect from (11121 the scaling behaviour [34, 35] 

M TRM (i, s) = ho [ du R(t, u) = r s- a f M {t/s) + r lS - x ^ z g M (t/s) , g M {y) ^ V~ Xr/z (6) 
Jo 

such that the asymptotic behaviour fmiv) ~ y~ Xl< / z holds true for y — > oo. In simple magnets, 
the leading correction term included in JHJ) often is quite sizeable and must be subtracted off 
before a reliable estimate of the scaling function /m(j/) can be obtained. 

By analogy with the simple Ising model, we determined the non-universal constants ro and 
T\ by fixing y = t/s. Then a fit was made from the plot of Mtrm(2/s, s) over against s. The 
results collected in table [I] have been obtained by averaging over at least three different values 
of y. As a first example, we show in figure^ the thermoremanent magnetization Mtrm(2/s, s) 
for e = 0.5 and T = 1. It turns out that in this case already the raw data scale very well which 
means that for the chosen values of the parameters, the finite-time corrections to the leading 
scaling Mtrm(*>s) ~ s_a /j\/(^/s) ar e negligible. We find the ageing exponent a — 0.398(5). 
In phase-ordering kinetics, the dynamics of systems with short-ranged spatial correlators is 
determined by the motion of the domain walls only [2]. Then simple dimensional analysis 
yields the relation a = 1/z. Comparing this with the values of a for which a data collapse is 
achieved, we find indeed a complete agreement with 

a = a{T,e) = ziT.e)- 1 (7) 

where z = z(T, e) is given by eq. ©, see the last column in table|l| We see that this agreement 
holds to within our numerical precision for all the values of T and e we have considered. In 
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Fig. 2 - (a) Thermoremanent magnetization Mtrm{us, s) as a function of s for e = 0.5 and T — 1.0 
and different values of y. The slopes of the fitted straight lines yield the average value 0.398(5) 
for the exponent a. (b) Time-dependent behaviour of the autocorrelator C(t, s = 0) and of the 
thermoremanent magnetization Mtrm^.s = 800) for e = 1 and T = 0.4. The slopes of the straight 
parts are different for large values of t, indicating that the exponents Xr/z and Xc/z are different. 
For illustrative reasons Mtrm has been multiplied by 10. 



this way, we confirm the earlier conclusion of Rieger et al. [30,31] on the scaling of the domain 
size L{t). We also recall that for the simple Ising model, one has a = 1/z = 1/2. ( 2 ) 

We point out that in contrast to the autocorrelation function, the integrated response 
displays a completely conventional dynamical scaling according to eqs. iHI2[) . 

Having understood how to achieve dynamical scaling for the linear response, we next 
look at the asymptotic behaviour /m(j/) ~ y~ Xli l z for y — > oo. We collect the results for 
the autoresponse exponent Xr/z in tableland now compare with the values of Xc/z we 
have determined before. It can be clearly seen that the autocorrelation and autoreponse 
exponents are different for the model at hand. We further illustrate this in figure^ where 
the autocorrelation function C(t, 0) ~ t~ x °/ z for t large enough and the thermoremanent 
magnetization Mthm(*j 800) ~ £--W z are compared for e — 1 and T = 0.4 in a log-log plot. 
Clearly the slopes for t large and hence the two exponents are different. This behaviour sets 
the kinetics of the bond-disordered Ising model apart from what is known for the phase- 
ordering of all simple magnets with short-ranged initial correlations. We did not anticipate 
this finding but recall that the random-phase sine-Gordon model also shows Ac ^ Xr slightly 
below its critical temperature [12]. 

After these preparations, we are ready to study in more detail the form of the scaling func- 
tion of the integrated linear response. It is straightforward to integrate the prediction 
(|2J of local scale-invariance with a = a' . We have (2-F1 is a hypergeometric function) 

fu{y) = y- XR,Z 2Fx ( 1 + a, ^ - a; ^ _ a + i ; I) ( 8 ) 



( 2 )It was claimed that because of dangerous irrelevant variables, the above-mentioned simple picture might not 
be valid in phase-ordering, leading to a ^ 1/z. For example in the 2D Ising model a = 1/4 was proposed [36]. 
Detailed studies of the simple 2D Ising model reconfirmed a = 1/z = 1/2 [18,20,34,37,38], in agreement with 
the standard picture and contra the claim of [36] . The support of eq. J7J in the diluted Ising model provides 
further evidence against the claim raised in [36]. 
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ln(t/s) 

Fig. 3 - Comparison of the numerically determined scaling functions /m(£/ s ) with the predictions 
coming from local scale-invariance (full lines), (a) e = 0.5 and T = 1, (b) e = 1 and T = 0.6, (c) 
e = 2 and T = 1. In the first case corrections to scaling are negligibly small. Error bars are smaller 
than the sizes of the symbols. 

In this expression, the values of the exponents and also the normalization are already fixed, 
see table H] We can now compare the data for the scaling function with the prediction © 
of LSI. In figure [3] examples for this comparison are shown for three choices of e and T. In 
the first case we have already seen that finite-time corrections to scaling were negligibly small 
and we now see from figure that the form of the scaling function is perfectly described by 
LSI. In the two other cases, the finite-time corrections are notable and actually quite sizeable 
for e = 2. Subtracting them, a very good scaling behaviour is found. Again, the curve © 
matches the data very well. This agreement with LSI could not hold if J5J or Q) were invalid. 

It is remarkable, in view of the very large values of z, how well the theory of local scale- 
invariance describes the data. From a field-theoretical perspective, it was advanced [4] that a 
symmetry principle such as local scale-invariance could at best hold at the mean-field level, 
where z = 2. The present example, using a non-integrable spin system with values of the 
dynamical exponent which can become very large, rather suggests that the idea of extending 
dynamical scaling to a more local form should indeed be capable of faithfully reproducing at 
least the linear responses of physically quite distinct systems. At present, it appears that the 
main requirements on the models to be studied are an algebraic growth law L(t) ~ t x / z of the 
linear size of the correlated clusters and a vanishing initial magnetization. The present study 
closes a gap between the simple, non-disordered magnets and on the other side critical spin 
glasses [14] where local scale-invariance has been successfully tested before. 

In summary, we have studied the ageing behaviour of a bond-disordered two-dimensional 
Ising model quenched to temperatures far below its critical temperature T c (e). While, the 
linear response as studied by the thermoremanent magnetization was seen to be completely 
compatible with standard dynamical scaling, large finite-time corrections to standard dy- 
namical scaling were seen in the two-time autocorrelation function. Our data appear to be 
compatible with the standard relation a = 1/z of the ageing exponent a with the dynamical 
exponent z(T,e) given by eq. ©■ We also encountered two surprises: first, assuming Q, 
the effective ageing exponent b was found to be non-vanishing (and negative) and second, 
the autocorrelation and autoresponse exponents are different, Ac 7^ Ar. In spite of the large 



Malte Henkel and Michel Pleimling: Disordered magnets & local scale-invariance 7 



range of values of z, the form of the scaling function of the thermoremanent magnetization is 
described to within the numerical accuracy by the theory of local scale-invariance. 

* * * 

MH thanks the Dipartamcnto di Fisica/INFN di Firenze for warm hospitality. This work 
was supported by Procope, by CINES Montpcllier (projet pmn2095), and by the Deutsche 
Forschungsgemeinschaft (grant no. PL 323/2). 

REFERENCES 

[1] L.C.E. Struik, Physical ageing in amorphous polymers and other materials, Elsevier (Amster- 
dam 1978). 

[2] A.J. Bray, Adv. Phys. 43, 357 (1994). 

[3] C. Godreche and J.-M. Luck, J. Phys. Cond. Matt. 14, 1589 (2002). 

[4] P. Calabrese and A. Gambassi, J. Phys. A38, R133 (2005). 

[5] M. Henkel, M. Pleimling and R. Sanctuary (eds), Ageing and the glass transition, Springer 
(Heidelberg 2006). 

[6] D.S. Fisher and D.A. Huse, Phys. Rev. B38, 373 (1988). 

[7] D.A. Huse, Phys. Rev. B40, 304 (1989). 

[8] A. Picone and M. Henkel, J. Phys. A35, 5575 (2002). 

[9] H.K. Janssen, B. Schaub and B. Schmittmann, Z. Phys. B73, 539 (1989). 

[10] A. Picone and M. Henkel, Nucl. Phys. B688, 217 (2004). 

[11] G.F. Mazenko, Phys. Rev. E69, 016114 (2004). 

[12] G. Schehr and P. Le Doussal, Phys. Rev. E68, 046101 (2003). 

[13] M. Henkel, Nucl. Phys. B641, 405 (2002). 

[14] M. Henkel and M. Pleimling, J. Phys. Cond. Matt. 17, S1899 (2005). 

[15] M. Henkel, T. Enss and M. Pleimling, J. Phys. A at press cond-mat/0605211 

[16] A.D. Rutenberg and A.J. Bray, Phys. Rev. E51, 5499 (1995). 

[17] C. Godreche and J.-M. Luck. J. Phys. A33, 9141 (2000). 

[18] M. Henkel and M. Pleimling, Phys. Rev. E68, 065101(R) (2003). 

[19] S. Abriet and D. Karevski, Eur. Phys. J. B37, 47 (2004); B41, 79 (2004). 

[20] E. Lorenz and W. Janke, to be published (2006); W. Janke in [5]. 

[21] M. Pleimling, Phys. Rev. B 70, 104401 (2004). 

[22] F. Baumann and M. Pleimling, J. Phys. A: Math. Gen. 39, 1981 (2006). 

[23] T. Enss, M. Henkel, A. Picone and U. Schollwock, J. Phys. A37, 10479 (2004). 

[24] J.R. Ramasco, M. Henkel, MA. Santos and CA. da Silva Santos, J. Phys. A37, 10497 (2004). 

[25] F. Baum ann, M. Henkel, M. Pleimling and J. Richert, J. Phys. A38, 6623 (2005). 

[26] G. Odor, cond-mat/0606724 

[27] H. Hinrichsen, J. Stat. Mech. Theor. Exp. L06001 (2006). 

[28] F. Baumann, S. Stoimenov and M. Henkel, J. Phys. A39, 4095 (2006). 

[29] C. Godreche, in [5] cond-mat/0604276 

[30] R. Paul, S. Puri and H. Rieger, Europhys. Lett. 68, 881 (2004); Phys. Rev. E71, 061109 (2005). 

[31] H. Rieger, G. Schehr and R. Paul, Prog. Theor. Phys. Suppl. 157, 111 (2005). 

[32] R. Paul, G. Schehr and H. Rieger, to be published (2006) . 

[33] A. Barrat, Phys. Rev. E57, 3629 (1998). 

[34] M. Henkel, M. Paefiens and M. Pleimling, Europhys. Lett. 62, 644 (2003). 

[35] M. Henkel, M. Paefiens and M. Pleimling, Phys. Rev. E69, 056109 (2004). 

[36] F. Corberi, E. Lippiello and M. Zannetti, Eur. Phys. J. B24, 359 (2001); Phys. Rev. E68, 

046131 (2003); Phys. Rev. E72, 028103 (2005). 

[37] C. Chatelain, J. Phys. A36, 10739 (2003). 

[38] M. Henkel and M. Pleimling, Phys. Rev. E72, 028104 (2005). 



